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Q\ , Abstract 
On 

We propose a parameterization for contributions of an infinite set of diagrams (bubble chains) into 

<-> 

physical observables represented as integrals of running coupling constant over the finite region in 



(N 



> 



momentum space. The perturbation theory part of contributions is rendered well-defined by the 
principal value prescrition for treating the Landau pole while the rest is connected to the gluon 
condensate in case of observables allowing the OPE analysis. The hypothesis of universality, i.e. the 



£S| ■ use of our parameterization for non-OPE cases, is discussed. 
(N 

in 

On 

q^. 1. Introduction. 

fH | Perturbation theory (PT) in the running coupling constant of strong interaction is fully understood 

and well developed technically being an excellent tool for exploring hard processes in QCD. However 
the numerical value of the effective parameter of the expansion is quite large [jjj that forces one to 
compute several first terms of perturbative expansion to get a reasonable precision. Experimental 
data are permanently getting better while new higher order terms of PT expansions appear rarely 
because of computational difficulties. This is especially discouraging because the standard model fits 
existing data well (2j and any sensible deflection will be noticeable at a high level of accuracy only. 

Aiming to precise comparison with experiment one should keep in mind that PT is asymptotic and 
cannot provide unlimited accuracy at fixed a s . In practice there is no indication yet on any asymptotic 
character of known series for physical observables, the main limitation being the technical complexity 
of getting new terms of an expansion. Few examples are known with several terms of a s expansion. 
For e + e~ annihilation M and the r semileptonic width 0] the second subsequent correction of order 
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is comparable with experimental precision in the MS scheme of subtraction. Terms of higher orders 
will hardly be available ever: perturbation theory seems to be saturated |J. 

And this is good because the computation of higher order terms is not worth it - being asymptotic 
the PT is also known to be incomplete. Theoretically recognized source of nonperturbative effects is 
the existence of instantons. Phenomenologically PT looks unsatisfactory as well - resonances can not 
be detected at any finite order in the coupling constant. 

Thus, stuck to QCD as a fundamental model it is reasonable to go beyond PT from inside by 
summing particular infinite sets of diagrams or specific contributions. For simple cases this has been 
done on several occasions. The necessity to sum infinite series appears, for instance, when an analytic 
continuation from Euclidean to Minkowskian region is performed H [7j. Some attention has been also 
drawn to fermionic bubbles summation in attempt to go beyond PT || (as a recent review, see ||). 

2. Fermionic bubble chains. 

The motivation for choosing this particular set of diagrams is threefold: 

i) in 1/Nf QED the corrections to photon propagator form a dominant set of graphs in Nf — > oo 
limit 0, 

ii) numerically the (3q expansion (a QCD counterpart of 1/Nf expansion) works well in MS sub- 
traction scheme for all known cases f\M , 

iii) they constitute a well defined set of graphs that can be summed. 

There are different and quite serious objections to above motivation but the third point is decisive; 
one wants to find diagrams that can be summed in a constructive way. 

Trying to keep this approach alive one could put forward some other foundations for fermionic 



bubbles chains at least as possible indicator of leading diagrams. Thus, skeleton expansions [12] or 
Schwinger-Dyson equations ][D| could lead in simple cases to the very same set of leading graphs 
in specific approximations. Early efforts to solve the confinement problem using a modified gluon 



propagator [14] lie in the same bed. 
3. Integral at small momenta. 

Whatever the underlying motivation is there are observables that can be represented on a general 
ground as an integral of running coupling including an infrared region. Examples have been discussed 
in the literature, we mention three of them e + e~ annihilation [||, pole mass of heavy quark [|l5|] and 
event shape [[u], [l7], |l8|| . Being infrared safe these observables can be formally written down as integrals 
over the running coupling that encounters the Landau pole present in the PT expression for a s (t) at 
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small t. Therefore their expansion in a s (Q), where Q is a large scale involved in the process, generates 
a factorial growth of PT coefficients. The emerging series is sign-definite and cannot be summed with 
Borel technique. Thus for a number of observables the representation of the form 

Q2 

a s (t)u>(t)dt (1) 



o 



can be formally obtained. The corresponding PT expression reads 

o f Q2 

a s (Q 2 ) / u(t)dt 



o 

and the difference between these two formulae reflects the result of "bubble summation improvement" . 
The weight function u(t) depends on kinematical invariants that are of order Q 2 . The problem reduces 
to a proper definition of the integral (|l]) that diverges due to the Landau pole if the perturbative 
running coupling is used in the infrared region of integration. 

Note that high orders of perturbation theory are normally analyzed using the Borel transform in 
the coupling constant. Then nonsummable pieces reveal themselves as poles in the complex plane of a 
Borel parameter. The growth of coefficients in fact can be much slower than factorial - it is enough for 
the series to be sign-definite to be nonsummable and to produce a singularity. Such singularities are 
called renormalons with specification from what region they come - IR or UV. So, factorial divergences 
are called renormalons though singularities can be connected with other types of large n behavior as 
well. In practice the simple way to generate coefficients at any n is to use the running coupling and 
the presence of Landau ghost is equivalent to existence of renormalons as factorial divergences of PT 
series. One could say that renormalization group sums logs while the proper definition of the behavior 
of the coupling constant at small momenta sums bubble chains. 

Qualitatively such chains are expected to produce nonperturbative pieces because a s (Q 2 ) and a s (t) 
are not connected perturbatively for t — ► and new terms are plausible, i.e. a new parameterization, 
different from purely perturbative, can enter the game. 

Dealing with factorial divergence in quantum field theory (rewritten as singularities in the Borel 
plane) at present state of art means to define how to treat the Landau ghost. Nevertheless in the 



literature the interpretation of the pole in the Borel plane is widely used IS]. For instance, the 



principal value prescription is very popular [2C]. In fact, the principal value prescription in the Borel 



plane is equivalent to the principal value prescription in momentum space at least at the level of one 
loop. Namely, a Borel represented observable of the form 



PV 



o 



1-t 



after the change of the variable t — > alog(Q 2 /t) goes to 

rQ 2 ( /-A 2 -e rQ 2 \ dt 

PV / a(t)dt = Um e ^o / + / -j — 7TT 

where we defined a(t) = (3oa s (t) = \jlogitj A 2 ). The last prescription is purely perturbative however. 

Namely, consider the integral with u(t) = 1, other functions do not change things qualitatively. Then 

rQ 2 rQ 2 rQ 2 SB. , x n 

PV / a(t)dt = / 2 a(t)dt = a(Q 2 ) / ^ ^ [oi{Q 2 )log{Q 2 /t)j dt (2) 

JO JQl JQl n=Q 

and 

rQl 

PV a(t)dt = 0, Ql = A 2 *1.45... 

The position of the zero Ql (Ql = A 2 * 1.45...) of the function 

rQ 2 

UQ 2 ) = PV / a(t)uj(t)dt 
Jo 

depends on the weight function uj(t) but the zero itself does exist for any function of the same sign at 
small t. The series @ is a convergent PT series. Indeed, inside the integration region we have 

\a{Q 2 )log{t/Q 2 )\ < \ a {Q 2 )logC^)\ = 1 - < 1. 

and the integrand converges homogeneously that allows one to integrate it getting again the convergent 
series. It is clear what happened - the most interesting region has been thrown away completely by 
choosing the PV prescription. 

Still this is a definition of the perturbation series. The coefficients of a s (Q) are not just numbers 
but are functions of Q 2 and contain powers of A/Q as well. In fact any cut will introduce another scale 
that renders the coefficients to become some functions of the ratio of those scales, in spirit though 
the expansion remains purely perturbative. In this sense the cut with principal value prescription 
is minimal because it does not introduce any new scales. The integral in eq. (||) can be computed 
explicitly 

I (log(Q 2 /t)) dt = Q 2 T(n + l,log(Q 2 /Q 2 )) 

jQo 

where T(n,z) is an incomplete T-function 

r z roo 

r(n + l,z)= e~H n dt = n\- / e~H n dt (3) 

Jo Jz 

that reduces to elementary functions for our particular case though. The last term in eq. (0) behaves 
asymptotically at large z as exp(-z), or exp(—l/a(Q*), where Q* = Q/y/lAh. Thus, any coefficient 
function of the new perturbative expansion (||) contains "nonperturbative terms" (see, also pl| ). 
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Numerically this recipe is valid meaning that two expressions 

r-Q 2 



PV / a(t)uj(t)dt and a(Q 2 ) / u{t)dt 
Jo Jo 



do differ. As we have shown however this difference can be accounted for perturbatively. We find 

PV f n Q2 a(t)u(t)dt 

ifi — -V w = 1.39, 1.27, 1.23 

a(Q 2 )/ Q u(t)dt 

for w(t) = 1 and Q 2 /A 2 = 100, 500, 1000. The convergence of the series looks like (for Q 2 /A 2 = 100) 

1.39 = 1 + 0.201 + 0.075 + 0.038 + 0.015 + . . . 

or in terms of a (100) 

1 + 0.92a(100) + 1.59a(100) 2 + 3.67a(100) 3 + 10a(100) 4 + 30a(100) 5 + . . . 

instead of factorial growth 

1 + l!a(100) + 2!a(100) 2 + 3!a(100) 3 + 4!a(100) 4 + 5!a(100) 5 + . . . 

The convergence is very slow. To reach a reasonable accuracy one needs almost as many terms as one 
could keep for the initial asymptotic series to get the best approximation. In this sense it imitates the 
asymptotic series very closely. 

The convergence can be essentially improved by choosing some other expansion parameter that 
reduces to the change of the scale a(100) — > a(100/£), £ > 1. The BLM choice [22] corresponds 



to vanishing of the first order correction; other optimization criteria lead to their own choice of the 
scale. On the whole, however, all these choices are perturbative in the sense that they give convergent 
series in the coupling normalized at some high scale. In practice convergence still can be slow and 
one is forced to use the integral formula but in principle this is a possible minimal way to define the 
perturbative series. Thus, the PV prescription does not create nonperturbative terms though it allows 
one to sum up some PT corrections in a closed form. Do the real nonperturbative terms exist in the 
chain? 

Look at the running coupling more carefully. The PV prescription suffers of being nonpositively 



defined at small momenta that can contradict some general properties of quantum field theory [23]. 
For instance, in the case of current correlators the spectral density must be positive. 
4. Extrapolation to infrared region. 
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We give several motivations and models for behavior of the coupling constant at small momenta 



keeping positivity p4\. 



First we use the freedom of choosing the renormalization scheme and coupling definition eventu- 
ally. Consider the e + e~ annihilation. Corresponding D(Q 2 ) function is known up to the order. 
Redefining a s in a scheme without higher order corrections we find 

D(Q 2 ) = 1 + ^1. 

7T 

Evolution of the new charge d s (Q 2 ) is governed by a new (3 function $(d s ) that has an IR fixed point 
and the effective charge is frozen at small momenta. The same conclusion has been recently obtained in 



[25] after analyzing the e + e~ annihilation cross-section R(s) within the principle-of-minimal-sensitivity 
approach [p6f |. Even without this reference one can choose a special (3 function providing a smooth 
infrared behavior for the coupling. For instance, 

a 2 

/3( a ) = "TT 2> K > ( 4 ) 

1 + kct 

where f3 as (a) = —a 2 + . . . p4[ . Note, that this extrapolation is based on PT consideration only and 
on the formal use of the summed form of PT series in a region where it is not supposed to be valid. 
One can think of eq. (||) as of a Pade approximation for the f3 function in a particular scheme. 

Some explicitly nonperturbative extrapolations exist also. The simple pattern is provided by the 
minimal analytic continuation 

a eff (Q 2 : 1 1 



that is regular everywhere outside a cut in the complex Q 2 plane but looks nonperturbatively in terms 
of an asymptotic charge a(Q 2 ) = (in 

a e ff(Q 2 ) = a (Q 2 ) - = a(Q 2 ) - e"^ + . . . 

and makes an extraction of "purely" nonperturbative terms untransparent. The corresponding (3 
function 

m = -e + - 1)- 2 

also has an explicit nonperturbative term unlike eq. (Q). 

Thus, there are some arguments in favor of smooth behavior of the coupling constant in the IR 
region keeping positivity. 

5. Nonperturbative terms and parameterization. 



Going back to our main problem of finding a difference between 

[Q 2 _ fQ 2 

/ a s (t)u(t)dt and a s {Q 2 ) I u(t)dt 
Jo Jo 

with a newly denned positive a s (t) one sees that new contributions should be added to the PV pre- 
scription - nonperturbative terms have appeared. It would be convenient to have a clear distinguish 
between PT and power corrections that are hidden in the exact formula (|l]). The simple parameteriza- 
tion could be also useful for comparison between different channels. We only need our function defined 
under the integration sign so we can account for the nonperturbative terms by localized distributions 
living at t = A 2 . We write p4| 

a e //(t) = a{t)\ pv + AA 2 5{t - A 2 ) + BA 4 5'(t - A 2 ) + . . . (5) 

Parameters A,B,... can be found for any particular model of extrapolation of running coupling 
constant into the infrared region. An important observation is that for a wide class of extrapolations 
an additional constraint A = can be satisfied. It is not necessary but looks plausible because in 
simple cases such a condition corresponds to the absence of gauge invariant operators with the mass 
dimension two. In other words we can extend our running coupling without the large distortion of 
PT. Several models belonging to the above class are: 
1. No interaction at small momenta 

a e ff(z) = a{z)Q{z - aA 2 ). 

System of equations for determining the cutoff a and parameters A and B is 

li(a)+A = 0, li(a 2 )-B = 

where li{a) is a special function 




with the PV prescription for the pole at real positive a > 1. A solution (with a constraint A = 0) is 
a = 1.45, B = «(2.1) = 1.19. 
2. Freezing (e.g. fill ) 

a eff (s) = a(aA 2 )0(aA 2 - s) + a(s)@(s - aA 2 ). 

The system of equations is 

* =U(a), -^- = U(a 2 )-B. 
Ln(a) 2ln{a) 
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with the solution a = 3.85, B = 2.6. 

3. Minimal subtraction 

aeff{s) = { MjfW) ~ VA^T )0(aA2 " S) + a(s)0(s " aA2) - (6) 

Again 

li(a) — ln(a — 1) = li(a), li(a 2 ) — a — ln(a — 1) = li(a 2 ) — B 

and a = 2, B = 2. 

4. Models generated by the @ function of eq. (H) 

The solution depends on the parameter k entering the expression for function (Q). For k ~ 2 
the solution is close to one of models (2) and (3) but for different k there may be no solution at all or 
one with a nonzero A. 

Thus, for reasonable extrapolations the solution subjected to the condition A = does exist and 
is stable enough, i.e. B does not change much. As for the model (3) one could allow the formula 
@ be valid for any s without any cutoff a. Then such an extrapolation, though quite legal, would 
be nonminimal in our sense and distort the PT strongly. It would introduce 1/Q 2 terms [27| in cases 
where they seem to be forbidden by operator product expansion. But in fact they were just correspond 
to a definition of the perturbative series. 

Now one can normalize our parameterization in a particular case to express the nonperturba- 
tive parameter B through some known quantities and to predict new ones. The place to turn for 



normalization is e + e annihilation where the operator expansion is known. We find [24] 

(—G 2 ) = ^-Bk\ (7) 

For (a s G 2 ) = (0.440 MeV) 4 the relation between the gluon condensate and A for this particular model 
is 0.330 MeV ~ 1.52B 1 / 4 A. For models (2) and (3) we find a constraint 0.330 MeV = (1.87±0.18)A M g. 
Taken literary it gives Aj^g = 180 db 20 MeV that is in a reasonable agreement with the present data. 
In general, our result means that the numerical values of the gluon condensate and of the parameter 
Ajyjg are compatible with each other for smooth continuation into the infrared region (like models 
(2,3)). For the model (1) numerics will be slightly different and in case of the model (4) it depends 
on k ranging from a bad solution through the standard values to some bad ones again. 

The generalization to two loop approximation for the asymptotic charge is straightforward. The 
Landau pole still exists (though its location is changed a bit) and the above machinery can be applied. 



8 



6. Universality. 



Here we assume the hypothesis of universality [18, 28, 29 1 to predict power corrections to a number 



of observables. We discuss the pole mass and the event shapes. In the framework of resummation of 
the perturbative corrections by a particular method of summation, the uncertainties, or the limit of 
accuracy, for physical observables are normally discussed. Our approach, connecting different channels, 
allows one to express one observable through another using the parameterization through the gluon 
condensate as an intermediate step and therefore gives the absolute magnitude of corrections. 

First example is the pole mass of a heavy quark. Being well defined in the PT framework it allows 
the representation of the form (||) by inserting the bubble chain into the one loop mass operator. The 
difference between the pole mass and the running mass entering the renormalized QCD Lagrangian is 



given by [15 1 



8vr f d 3 k a s (k) 2 f» a s (t)dt 



Using our parameterization (@) we get the result 



B A AB 
mp — run = A = 



at /3q = b/Air which should be compared with the uncertainty found in [15 



m P -m Q = —A. 

The value of the parameter B is determined by the gluon condensate and is about two from (|^). 
Second example is the mean value of the thrust parameter T that has been computed up to the 



next-to-leading order in a s [ 30 1 



(T) = 1 — 0.355a s (Q 2 )(l + 9.56^-^). 



7T 



Corresponding uncertainties due to bubble chain summation and estimates of generated power correc- 
tions have been considered in Rig , |l7], [l£|. We use the representation of the form (jlj) given explicitly 



in |7| 



16 fQ „ . 8 [Q 2 a s (t)dt 

Our parameterization leads to the result 

Note, that in this particular example formula @ coincides with the previous case (||) up to numerical 
factors. 



The dependence on the scheme and specification of A come from the analyses of the perturbative 
part of the expansion. If the particular scheme is chosen for PT expansion then the corresponding A 
appears in @. We are going to elaborate on this point elsewhere. 

7. Conclusion. 

We propose to parameterize the power corrections in a well-defined way. The universality hypoth- 
esis (physical justification of which we don't discuss here) allows one to connect these corrections in 
different channels. The PT structure is reflected in the choice of the parameter A: as soon as the PT 
scheme is fixed the corresponding parameter A appears in power corrections in our approximation. 
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a part of this paper was written and interest in the work. AAP thanks K.G.Chetyrkin, M.Jezabek 
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